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Abstract 

Motivated by the recent surprising discovery of two new meson states with cs quark 
content but unexpectedly low masses and narrow total decay widths, we work out, in a 
nonrelativistic potential-model approach developed already some two decades ago, the 
predictions for the energy levels of the corresponding charm-anti-strange quark bound 
states. In spite of the fact that this simple quark model reproduces the mass spectrum 
of the previously observed hadrons remarkably well, we are led to the conclusion that, 
without considerable modifications, both the new states do not fit into this framework. 
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1 Introduction 

Recently, evidence for the existence of two new mesonic states with charm and strange 
quark content has been presented by three experiments [1-13]. Table 1 summarizes the 
experimental measurements of the masses of the two new states as well as our weighted 
averages of the meson masses. The uncertainties of the averages have been obtained by 
adding in quadrature the statistical and systematic errors reported by the experiments. 

At first, the BaBar Collaboration [1], running the BaBar detector at the PEP-II 
asymmetric-energy e+e~ storage ring, announced the observation of a charmed strange 
meson called D*j{2317)^. This new state manifested as narrow resonance with a mass 
ofM(L);j(2317)±) = (2316.8±0.4±3.0) MeVinthe DfTr'^ invariant-mass distribution. 

The CLEO Collaboration [2,3, 12] provided the confirmation of the existence of the 
D*j{2317)^ meson from data collected with the CLEO II detector in symmetric e+e~ 
collisions at the Cornell Electron Storage Ring. Furthermore, the CLEO Collaboration 
[2,3,12] reported the observation of a charmed, strange meson called Dsj{2463)^. This 
new state manifested itself as narrow resonance in the D*^7r^ invariant-mass spectrum. 
Interestingly, the CLEO experiment [3,12] found, within errors, identical results for the 
mass differences (with statistical and systematic errors) between these two D*j{2317)^ 
and Dsj{24:63)'^ states and the well-estabhshed [14] Df and D*^ mesons, respectively: 

M{Dlj{2317)^) - M{Df) = (350.0 ± 1.2 ± 1.0) McV , 
M(D,j(2463)^) - M{Df) - (351.2 ± 1.7 ± 1.0) MeV . 

With the most recent Particle Data Group averages [14] M{Df) = (1968.1±0.5) MeV, 
for the mass of the Df meson, and M(L)*=^) = (2111. 9±0. 7) MeV, for the mass of the 
D*^ meson, the mass differences translate into the resonance masses listed in Table 1. 
In return the BaBar experiment confirmed the existence of the Dsj(2463)^ meson [7] . 

Finally, both discoveries have been confirmed by the Belle Collaboration, operating 
the Belle detector at the KEKB asymmetric-energy e"'"e~ collider; the Belle experiment 
identified the new L'*;(2317)=^ and Dsj{24Q3)^ states in the exclusive B meson decays 
B D+D* j{2317) and B D+Dsj{2A63) [5,11], and continuum e+e" annihilations 
[6]. 



Table 1: Masses of the two new narrow charmed, strange meson states £'*j(2317)^ and 
Dsj{2A63)^ with both statistical and systematic errors, measured by the experiments 
BaBar [1,7], CLEO [2,3,8,12], and Belle [5,6,9,11], as well as their weighted averages 



Experiment 


Process 


M(D*j(2317)±) [MeV] 


M(L>,j(2463)±) [MeV] 


BaBar [1,7] 


inclusive c+c^ 


2316.8 ±0.4 ±3.0 


2457.0 ±1.4 ±3.0 


CLEO [3,8,12] 


inclusive c+c^ 


2318.1 ± 1.2 ± 1.0 


2463.1 ± 1.7 ± 1.0 


Belle [5, 11] 


B decays 


2319.8 ±2.1 ±2.0 


2459.2 ±1.6 ±2.0 


Belle [6] 


inclusive e+e" 


2317.2 ±0.5 ±0.9 


2456.5 ±1.3 ±1.1 


Weighted Averages 


2317.6 ±0.8 


2459.0 ± 1.1 
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Until now, for the L'*j(2317)^ state, only the decay D*j{2317y D++7i^ and, for 
the D,j(2463)^ meson, the decays Dsj{2A63)+ 0*+ + %^^ and D,j(2463)+ ^ D+ + -f 
(as well as, of course, their charge-conjugate channels) have been observed [1-9,11,12]. 
For both states, the widths of the resonance peaks are consistent with the experimental 
mass resolution. This indicates that the observed resonance widths are consistent with 
those expected for states of zero intrinsic or natural width F. The CLEO Collaboration 
reports identical upper limits: F(D*j(2317)^) < 7 MeV for the D*j{2317)'^ meson and 
F(D,j(2463)±) < 7 MeV for the Dsj{2m)^ meson, at the 90% confidence level [3,12]. 

Concerning the quantum numbers of these two narrow resonances, all experimental 
findings, such as all the observed decay modes and angular distributions, are consistent 
with their interpretation as P-wave states with spin-parity assignment = 0^ for the 
D*j{2317)^ meson [1-4,6-9,12] and = 1+ for the L',j(2463)± meson [2-9,11-13]. 

The new states do not show the decay patterns favoured by theoretical expectation: 

• The mass of the L'*j(2317)^ meson is 40.8 MeV below the kinematical threshold 

2358.4 MeV for its isospin-conserving strong decay £)*j(2317)+ D° + K^, and 
49.4 McV below the kinematical threshold 2367.0 McV for its isospin-conserving 
strong decay D*j{2317)~^ D~^ + K^ (charge-conjugate decays are understood). 

• The mass of the Dgj(2A63)^ meson lies above the two kinematical thresholds for 

the decays L',j(2463)+ -> D°+K+ and D,j(2463)+ ^ D++K^; these processes, 
however, are in conflict with the conservation of angular momentum and parity. 

• The mass of the Dsj(2A63)^ meson is 41.5 MeV below the kinematical threshold 

2500.5 McV for its isospin-conserving decay D,j(2463)+ D*{2007f+K+, and 
48.7 McV below the kinematical threshold 2507.7 MeV for its isospin-conserving 
decay Dsj(2463)+ D*{2010)^+K^ (charge-conjugate decays are understood) . 

As a consequence of their low masses blocking their isospin-conserving strong decays, 
the mesons D*j{2317)^ and L>,j(2463)± decay via the channels D*j(2317)+ Df+n° 
and L'sj(2463)+ — > D*+-|-7r°, which both violate isospin conservation, and through the 
radiative decay mode Dsj{2463)~^ D^+'-f, as well as the respective charge- conjugate 
modes. This explains the observed narrowness of the total decay widths of these states. 

The experimentally measured masses of the two D*j(2317)^ and Z)sj(2463)^ states 
(summarized in Table 1) are significantly smaller than the corresponding predictions of 
most, but not all, QCD-inspired quark potential models, which regard these mesons as 
P-wave bound states (cs) of a charm quark c and a strange antiquark s. Concerning the 
masses of the = 0^ and the lower = 1+ (cs) states, typical sets of predictions are 
obtained, e.g., in the comprehensive potential-model analyses of Refs. [15-17], whereas 
results at the low-mass end of the spectrum of predictions can be found in Refs. [18, 19] . 
The failure of quark potential models to explain such low masses of the new resonances 
prompted a multitude of speculations about the genuine nature of these states [20-58]. 
Tentative interpretations view the mesons D*j{2317)^ and Z^sj(2463)'^ dominantly as 
(cs) bound states [20,22,24,26-28,31,36,38,46,49,51], (DK) molecules [21], four-quark 
bound states [23,29,33,53], or Dn atoms [25]. Clearly, in general each of these mesons is 
a linear combination of all Fock states which carry the appropriate quantum numbers. 

In this work, we examine whether the new Dg mesons can be accommodated as (cs) 
bound states in a phenomenologically appealing nonrelativistic quark potential model. 
After recalling in Sec. 2 the spectroscopic classification of heavy-light states, we sketch 
in Sec. 3 the adopted quark model, and present in Sec. 4 our — disappointing — findings. 
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2 Spectroscopy of heavy-quark light-quark mesons 

The spectroscopic classification of bound states (Qq) composed of a heavy quark Q and 
a comparatively light antiquark q is greatly facilitated by the following observation. In 
the limit of infinitely large mass of the heavy quark Q, the spin Sq of the heavy quark 
Q and the total angular momentum j of the light quark q, defined as the sum j = l+Sq 
of the orbital angular momentum I and the spin Sq of the light quark q, are separately 
conserved. Hence, in this case j provides a "good" quantum number for classification. 

The light-quark angular momentum j and the spin sq of the heavy quark combine 
to the total angular momentum J = j+Sq of the (Qq) bound state. The corresponding 
quantum number J fixes the spin of the resulting meson. The parity P of this meson is 
related to the quantum number i of the orbital angular momentum Z by P = (— l)^"*"^. 

Thus, in the heavy-quark limit the heavy-light bound states (Qq) may conveniently 
be classified in terms of the light-quark total-angular-momentum quantum number j: 

S-wave states have orbital angular momentum i — and, therefore, negative parity, 
P = —1. Spin and orbital angular momentum of the light quark q can couple only 
to J'^ = I • Combining the light-quark total angular momentum j = Sg with the 
spin of the heavy quark Q yields a spin-singlet state with spin-parity assignment 
= 0~ and a spin-triplet state with spin-parity assignment = 1~ . 

The pseudoscalar = 0^ state has been identified with the isosinglet Df meson 
(the former meson), with a mass of (1968. 1±0. 5) MeV and a well-established 
spin-parity assignment = 0~ [14]. The vector = 1~ state is assumed to be 
identical to the isosinglet D*^ meson, with a mass of (2111. 9±0. 7) MeV, natural 
spin-parity, and width and decay modes consistent with the assignment = 1~ 
[14]. 

P-wave states have orbital angular momentum i = 1 and, therefore, positive parity, 
P — +1. Spin and orbital angular momentum of the light quark can couple either 
to — 1^ or to — 1^. Combination with the spin of the heavy quark yields, 
for — 1^, two states with spin-parity assignment — 0"^ and — and, 
for — I"*", two states with spin-parity assignment — 1+ and — 2+. The 
two — l"*" states do not have definite charge-conjugation properties; therefore, 
they can undergo mixing. 

The vector = 1+ state belonging to the = |^ doublet (with possibly small 
admixtures of its vector = l"*" counterpart belonging to the = I"*" doublet) 
is, in general, assumed to be identical to the isosinglet Dsi{2536)^ meson, with a 
mass of (2535.35±0.34±0.5) MeV and its spin-parity assignment J^' = 1+ which 
is strongly favoured but still needs confirmation [14] . 

The spin-triplet tensor — 2+ state is identified with the isosinglet Dsj{2573)^ 
meson, with a mass of (2572.4 ± 1.5) MeV, natural spin-parity, and width and 
decay modes consistent with the assignment — 2+ [14]. 

Therefore, the present phenomenological status of the Ds system can be summarized in 
the following way: The two S-wave (ground) states are experimentally well established. 
Good candidates for two of the four P-wave states predicted by theory for bound states 
of two spin-| fermions have been observed. The open question is: May the D*j(2317)^ 
and Dsj{2A63)^ states be interpreted as one or the other of the two missing (cs) states? 
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3 Nonrelativistic quark— ant iquark potential model 

The present analysis of the charmed strange meson system is based on a nonrelativistic 
quark model improved by the lowest-order relativistic corrections. For recent reviews of 
the phenomenological description of hadrons as bound states of quarks within the more 
intuitive framework of (both nonrelativistic and semirelativistic) potential models see, 
for instance, Refs. [59,60] . Since we intend to describe bound states composed of both a 
heavy quark and a light quark, we have to make use of a framework which is capable to 
reproduce the properties not only of heavy quarkonia but also of light mesons [61-63]. 

The nonrelativistic interquark potential V{r) characterizing the phenomenological 
model employed in the present analysis is of a "funnel- like" Coulomb-plus-linear shape, 
improved by an intermediate exponential damping [61]. This quark model provides an 
excellent description of both the observed [14] meson and baryon mass spectra [62,63]. 
Apart from the strong finc-stnicture constant and the masses of the qiiarks, this 
model involves six constant parameters, a, b, c, d, Vq, Ri, entering in the potential V{r): 

V(r) = -- — -6exp — )+d + Vo for r < i?i , 
3 r \ cj 



V{r) = ar — bexp (^—-^ + Vq for r > it!i . 



Prom the smoothness of V{r) at r = the potential parameters Ri and d are related 
to the strong fine-structure constant ccg and the slope a of the linear term according to 



Starting from the Bethe-Salpeter formalism for the description of bound states within 
relativistic quantum field theory, the static potential V{r) in the Schrodinger equation 
may be derived from the interaction kernel entering in the Bethe-Salpeter equation by 
Fourier transformation. However, for a confining potential, such as the linear potential 
V{r) = ar, the corresponding Bethe-Salpeter kernel has to be a highly singular object. 
The necessary regularization of the kernel introduces a constant in the static potential. 
In order to compensate the (infrared) divergence of the interaction kernel, in the case of 
the linear potential V{r) — ar this (negative) constant Vq has to involve the slope a of 
the linear potential, and the Euler-Mascheroni constant 7 = 0.577 215 664 901 . . . [64] : 

Vo = -2Vaexp [-(7 - |)] ~ -2^0 • 

The above relations reduce the number of parameters in the potential to four: ctg, a, b, c. 
With respect to the Lorentz structure of the quark-ant iquark interaction, the potential 
V{r) is split into a Lorentz-scalar part Vs = ar and a Lorentz- vector part Vy = V — Vs- 
The generalized Breit -Fermi Hamiltonian contains all relativistic corrections to the 
(nonrelativistic) Schrodinger Hamiltonian up to and including the second order in the 
inverse masses of the bound-state constituents. This lowest-order relativistic correction 
discriminates between different Lorentz structures of the quark-ant iquark interaction. 
The spin-dependent relativistic corrections consist, for a Lorentz- vector spin structure, 
of a spin-orbit, a spin-spin, and a tensor term, and, for a Lorentz-scalar spin structure, 
of only a spin-orbit term. The spin-spin term, however, involves the Laplacian AVv(^) 
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of the vector potential Vy{r); for any Coulomb-like contribution to the vector potential 
Vv{r), i.e., a term in Vv(r) proportional to 1/r, this yields a highly singular expression: 

A- = -47r5(^)(r) . 
r 

Similarly, for a Coulomb- like contribution V{r) oc 1/r to the nonrelativistic interaction 
potential V{r) the spin-orbit and tensor terms introduce singularities of the form 
This makes the Breit-Fcrmi Hamiltonian an operator which is unbounded from below. 
However, for heavy bound-state constituents, all relativistic corrections may be treated 
as perturbations of the well-defined nonrelativistic Schrodinger Hamiltonian operator. 

The mass difference between the pseudoscalar tt and the vector p meson, generated 
by the spin-spin interaction, is of the same order of magnitude as the involved meson 
masses. This indicates that, for any description of mesons containing light quarks, the 
inclusion of at least the spin-spin term in the nonperturbative treatment is mandatory. 
Within the nonrelativistic framework, a conceivable solution is the ad-hoc replacement 
of the 5 function by a suitable smearing function /(r) which converges weakly towards 
the S function in the appropriate limit, i.e., the use of a representation of the 5 function. 
This procedure should be independent of the particular choice of the representation of 
the S function. The model of Ref. [62] replaces the S function by the smearing function 

f(r) = ^ exp (-r/rp) 
•'^^ 47rr2 r 

satisfying 

hm /(r) = 5^^\r) . 

The dependence of the smearing-range parameter introduced here, Tq, on the masses of 
the bound-state constituents is encoded in a purely phenomenological relationship [62] . 
The numerical values of the parameters of this model have been determined in Ref. [62] 
from a fit of the mesons known at that time, yielding, for the constituent quark masses, 
mu = ma = 336 MeV, = 575 MeV, = 1845 MeV, = 5235 McV, and, for the 
couplings in the potential, = 0.31, a = 0.15 GeV^ b = 0.956 GcV, c = 2.05 GcV"^ 

For this quark interaction, the mass eigenvalues and eigenstates of the Brcit-Fermi 
Hamiltonian are computed to high accuracy with the help of the numerical integration 
procedure for the solution of the Schrodinger equation developed in Ref. [65]. However, 
before applying the above quark model to the main target of the present investigation, 
viz., the system of charmed strange (Dg) mesons, we would like to confront this model's 
predictions for the entire mass spectrum of the experimentally well-established mesons 
with results obtained in some representative different approaches. For this comparison, 
we choose the relativized quark potential model of Godfrey, Isgur, and Kokoski [15,16], 
and a relativistic quark model for mesons developed over the last decade by some Bonn 
group within the Lorentz-covariant framework of the Bcthc-Salpeter equation [66-71]. 

In order to get some idea of the significance of our approach, that is, to estimate the 
accuracy of our predictions, and to compare them with the findings of different models, 
we introduce, as a measure of quality, a quantity Q, defined as the average of the square 
of the relative error of the theoretical result m/*^^ for the mass of the composite particle 
i — 1,2,3,... with respect to the corresponding experimentally measured value M'f^^^ : 




^(exp) 
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The smaller the numerical Q value, the higher is the accuracy of the adopted approach. 
For the numerical comparison of the models, we take into account their predictions for 
N — 2A energy levels corresponding to the well-estabhshed mesons tt, p(770), 0(1020), 
7r(1300), p(1450), K, X*(892), D, D*(2010), D„ Dl B, B\ B^, Vc{lS), W(1S), 
Vc{2S), V'(2S), ^(4040), ^'(4415), T(1S), T(2S), T(3S), and T(4S). Numerically we find 
for the quality measure Q in our approach Q = 0.019%, for the Godfrey-Isgur-Kokoski 
model [15,16] Q = 0.041%, which is roughly twice as large as our Q value, and from the 
results presented in Refs. [68,70] for the Bethe-Salpeter treatment Q — 0.28%, which is 
more than an order of magnitude larger than our Q value. Prom these numbers, we feel 
entitled to conclude that this quark potential model is definitely competitive (at least). 

4 Results, Discussion, and Conclusions 

Table 2 gives the predictions of the quark model sketched in Sec. 3 for the masses of the 
lowest-lying six (cs) bound states expected by the spectroscopic classification of Sec. 2. 
Notice that we still adopt the model parameter values determined in the fit of Ref . [62] . 
The formalism of the Breit-Fermi Hamiltonian is more closely related to L-S coupling. 
Consequently, the states are labelled here by the usual spectroscopic notation 
where n denotes the principal quantum number, S is the quantum number of the total 
spin S — Sq + Sq, and S, P, D, . . . indicate orbital angular momentum ^ = 0,1,2,.... 
There are two — 1+ states: the spin-singlet ^Pi state and the spin-triplet ^Pi state. 
Because of non-diagonal contributions of the rclativistic corrections in the Breit-Fermi 
Hamiltonian, the two = 1^ states do not constitute mass eigcnstatcs. The physical 
states associated to the lower (1) or the higher (h) mass eigenvalue, respectively, of the 
Breit-Fermi Hamiltonian are called D^g} and D^gl\ The latter states are obtained from 
the states ^Pi and ^Pi by an orthogonal transformation involving some mixing angle 9: 

= cos^l^Pi) -sin^l^Pi) , 
IdJ?) = sin^l^Pi) +cose|^Pi) . 

In our quark model we find for this ^Pi-^Pi mixing angle the numerical value 9 — 44.7°. 

Moreover, Table 2 compares the predictions of our approach with, on the one hand, 
the experimentally observed masses [14] of all previously known Ds mesons and, on the 
other hand, the results of the four (nonrelativistic) potential models [15-19] explicitly 
mentioned in the Introduction and of the Bonn Bethe-Salpeter model [66-71]. First of 
all, we note that all predictions of the Lorentz-covariant Bethe-Salpeter framework do 
not substantially differ from the potential-model ones; this indicates that the problems 
with the new states are not just a consequence of the nonrelativistic approximation. 
Our results show a satisfactory agreement with experiment and arc within the range of 
predictions of the other models. We are forced to conclude that the model of Ref. [62] , 
too, cannot explain the masses of the new Dg resonances; this discrepancy becomes still 
more distinct when inspecting the binding energy of the quark-antiquark bound states. 

We may get a clue to the correct identification of the heavier of the two new states, 
the Dsj{2463)'^ meson, by inspection of the center of gravity (COG) of the £ = 1 states. 
For P-wave states, the mass M(COG)p of the center of gravity is defined according to 



M(COG)p = ^ [M(^Po) + 3 M(3Pi) + 5 Mi^P^)] 
9 
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Table 2: Energy levels (in units of MeV) of the lowest S- and P-wave states, identified 
by the usual spectroscopic notation of the charmed, strange quark-antiquark 

{Dg meson) system predicted by the nonrelativistic quark model developed in Ref. [62]; 
the two mass eigenstates of the two-state = 1+ P-wave subsystem corresponding to 
lower and higher mass, formed by linear combinations of the spin-singlet ^Pi state and 
the spin-triplet ^Pi state, are denoted by Dg} and D^^^ , respectively. For comparison, 
the experimental values and some results of typical quark models [16-19] arc listed too. 
Since Refs. [16] and [19] round their results to 10 MeV zeros are added where necessary. 
The ^,1(2536)^ of mass (2535.35±0.6) MeV is one of the = 1+ states and D^^. 



State 


Experiment 


This Work 


Ref. [16] 


Ref. [18] 


Ref. [19] 


Ref. [17] 


Ref. [70] 


I'So 


1968.1 ± 0.5 
2111.9 ±0.7 


1963 
2099 


1980 
2130 


1968.8 
2110.5 


1940 
2130 


1965 
2113 


1969 
2116 


l=^Po 

1^P2 


2572.4 ± 1.5 


2446 
2515 
2527 
2561 


2484 
2550 
2560 
2590 


2387.8 
2521.2 
2536.5 
2573.1 


2380 
2510 
2520 
2580 


2487 

2535 
2605 
2581 


2464 
2506 
2506 
2552 



Within the framework of (nonrelativistic) interaction-potential models for two-particle 
bound states, a perturbative inclusion of the spin-spin interaction implies the equality 
of the P-wave center-of-gravity mass M(COG)p and the mass M(-^Pi) of the ^Pi state: 

M(COG)p = M(^Pi) (perturbatively) . 

In contrast to this, in a nonperturbative treatment of the spin-spin interaction by, e.g., 
some smearing of Dirac's 5 function — which should come closer to reality — the P-wave 
center-of-gravity mass M(COG)p will be larger than the mass M(^Pi) of the ■'^Pi state: 

M(COG)p > M(iPi) (nonperturbatively) . 

Let's neglect the mixing of the two = 1^ states necessary to get the physical states. 

• Identifying the ^Pi state with the L),,i(2536)^ of mass (2535.35±0.34±0.5) MeV, 
and the ^Pi state with the Dsj{2A63)^ meson of mass (2459. Oil. 1) MeV, entails 

M(COG)i,'"P^ = 2531.7 MeV > M(^Pi) = M(D,j(2463)^) = 2459.0 MeV , 
in accordance with the above general theoretical expectations for this inequality. 

• Identifying the ^Pi state with the Dsj(2463)^ meson of mass (2459.0±1.1) MeV, 
and the ^Pi state with the Dsi(2536)± of mass (2535.35±0.34±0.5) MeV, entails 

M(COG)J,'''P^ = 2506.3 MeV < M(^Pi) = M(D,i(2536)±) = 2535.35 MeV , 

in contradiction to the above general theoretical expectation for this inequality. 

This simple observation may be regarded as a hint that the new i^sj (2463)^ resonance 
is predominantly a ^Pi state. 
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